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Abstract 

Assuming that the lightest neutral component in an SU(2)l gauge multiplet is 
the main ingredient of dark matter in the universe, we calculate the elastic scattering 
cross section of the dark matter with nucleoli, which is an important quantity for the 
direct detection experiments. When the dark matter is a real scalar or a Majorana 
fermion which has only electroweak gauge interactions, the scattering with quarks 
and gluon are induced through one- and two-loop quantum processes, respectively, 
and both of them give rise to comparable contributions to the elastic scattering cross 
section. We evaluate all of the contributions at the leading order and find that there 
is an accidental cancellation among them. As a result, the spin-independent cross 
section is found to be 0(1(H 46 - 48 ) ) cm , which is far below the current experimental 
bounds. 



1 Introduction 



The existence of dark matter (DM) is one of the mysteries of the universe. Its energy 
density in the universe, which is about six times larger than that of baryon [1], cannot 
be explained in the standard model (SM) in particle physics. Weakly interacting massive 
particles (WIMPs) beyond the standard model are regarded as promising candidates for 
the DM. If they have TeV-scale masses, their relic abundance in the thermal history of 
universe may naturally account for the observed value. For the past years, a lot of efforts 
have been dedicated to the direct detection of WIMP DM, and their sensitivities have 
been improving. The XENON100 experiment [2], for example, has already started and 
announced its first result, which gives a stringent constraint on the spin-independent (SI) 
WIMP-nucleon elastic scattering cross section afj ( a^} < 2.4 x 10 -44 cm 2 for WIMPs with 
a mass 50 GeV [3]). Furthermore, ton-scale detectors for the direct detection experiments 
are now planned and expected to have significantly improved sensitivities. 

Introduction of new particles with masses TeV scale is one of the simplest extensions 
of the SM in order to explain the DM in the universe if they are an electrically neutral 
component in an SU(2) L gauge multiplet (n-tuplet) with the hypercharge Y. We call 
them electroweak-interacting massive particles (EW-IMPs) in this article. EW-IMPs are 
assumed to interact with quarks and leptons only via weak gauge interactions. Neutral 
Wino and Higgsino in the minimal supersymmetric standard model (MSSM) are examples 
of such particles when the mixing with other states is negligible and squarks and slep- 
tons are heavy enough. There is an earlier work |3] which studied the direct detection of 
Higgsino-like DM. They gave the cross section of Higgsino-like DM with nucleon by point- 
ing out large loop contribution to Yukawa couplings. The cross sections of Wino/Higgsino 
DM with nucleon were calculated in Ref. [5]. In their study, they took into account elec- 
troweak loop corrections and showed that the SI interaction does not vanish in the large 
DM mass limit, which is a distinctive feature of EW-IMP DM. There are several other ar- 
ticles which give theoretical prediction for the cross section of EW-IMP DM. The authors 
in Ref. [6] intensively studied EW-IMP DM (they refer to it as Minimal dark matter), 
and a similar calculation was also performed in Ref. [TJ. There is, however, inconsistency 
among Refs. [6J E] - In addition, interaction of gluon with EW-IMP, which is also one 
of the leading contributions in the EW-IMP DM-nucleon scattering, was neglected in the 
references. This point was first pointed out in Ref. [8]. 

In this paper, we give accurate prediction for EW-IMP-nucleon scattering in the direct 
detection experiments. Assuming that EW-IMPs are the main component of the DM in 
the universe, we calculate the cross section at the leading order of SM gauge couplings, 
including the gluon contribution. We provide the complete formulae for the EW-IMP 
DM-nucleon scattering cross section in a general form, and show that the SI cross section 
is below the current experimental bounds in all the cases we have studied. We also give 
numerical results for spin- dependent (SD) cross section. 

This paper is organized as follows: in Section [2] we explain the EW-IMP DM scenario. 
In Section [3] general formulae of the elastic scattering cross section of dark matter with 
nucleon are summarized. In Section H] we derive the cross section of the EW-IMP with 



1 



nucleoli and show the numerical results. Section is devoted to conclusion. 

2 EW-IMP Dark Matter 

As we described in the Introduction, EW-IMPs are an electrically neutral component in 
n-tuplet of SU(2)l with the hypercharge Y. In this section, we explain characteristics of 
EW-IMP DM assuming the EW-IMPs are fermionic. Fermionic EW-IMPs are popular, 
e.g., Wino or Higgsino in the MSSM. We will discuss scalar EW-IMPs in the end of this 
section. 

When the fermionic n-tuplet has only SU{2)l x U(1)y gauge interactions in the SM 
and all the components have a common mass at tree level, the charged components become 
heavier than the neutral one because of quantum loop corrections [B]. The typical value 
of the mass difference AM is 0(100) MeV. Thus, the neutral component of the multiplet 
could be a DM candidate in the universe if it is stable. It is found that fermionic EW- 
IMPs with n > 5 could be stable without imposing a certain symmetry to forbid its 
decay, such as R parity in the MSSM, since the gauge and Lorentz invariance prevent the 
DM candidate from decaying via the renormalizable interactions, and they turn out to 
be stable even if non-renormalizable effective operators are considered under the cut-off 
scale as large as the GUT-scale or Planck-scale [6j 9_|E 

The thermal relic abundance could also explain the observed DM energy density when 
the EW-IMP mass is over 1 TeV for n > 2. In Ref. [TO], the thermal relic abundance 
of Wino in the MSSM (which is EW-IMP with n = 3 and Y = 0) is evaluated^ The 
thermal relic abundances for EW-IMPs with n = 2 and 5 (Y — 0) are also evaluated [9]. 
According to those studies, the EW-IMP mass 1, 2.7, and 10 TeV is suggested in order 
to explain the observed DM abundance when n — 2, 3, and 5, respectively, and Y = 0. 

Now we explicitly show the Lagrangian of EW-IMP DM scenario for our calculation. 
In the case of the multiplet with the hypercharge Y — 0, the following Lagrangian is 
introduced to the SM : 

A£ = ±Z{ip-M Q )x, (1) 

where x is an SU{2)l n-tuplet fermion, and _D M = <9 M — igiYB^ — ig^W^T 0, [a = 1, 2, 3) 
is the gauge covariant derivative. Here g\ and g<i is U{1)y and SU{2)l gauge coupling 
constants, respectively, and T a is for the generators of SU{2)l gauge group. The tree- 
level mass of the multiplet is denoted as M . The neutral component of the n-tuplet is a 
Majorana fermion, while the other charged components are Dirac fermions. The neutral 
component has no interaction with gauge bosons by itself so that the elastic scattering of 
EW-IMP with nucleon is induced via loop diagrams. 

1 In this article we implicitly suppose a symmetry to prevent EW-IMPs from decaying for 2 < n < 5. 
2 In the evaluation of thermal relic abundance of EW-IMPs with n > 2, the Sommerfeld effect in the 
annihilation cross section of EW-IMP [TTJ should be included [TU] ■ 
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On the other hand, for Y ^ n-tuplet, the following term is added to the SM La- 
grangian, 

AC = - M )^ . (2) 

Here, is an SU(2)l n-tuplet Dirac fermion so that the fermion have an SU(2)l x 
U(1)y invariant mass term. Note that Dirac fermions as the DM in the universe have 
been severely constrained by the direct detection experiments because of its large elastic 
scattering cross section via coherent vector interaction with target nuclei. 

The situation changes if introducing the extra effective operators which give rise to 
a mass splitting 5m between the left- and right-handed components of the neutral Dirac 
fermion after the electroweak symmetry breaking. In this case, the neutral component is 
no longer mass eigenstate; the mass splitting decomposes the neutral Dirac fermion, 
into two Majorana fermions, x° and fp : 

^ = ^=(x° + iff). (3) 

Hereafter we treat x° as the lighter one, without loss of generality. Now the lighter 
component could be the DM candidate. It does not have any vector interaction by itself 
so that the DM-nucleon elastic scattering by Z-boson exchange is forbidden. In addition, if 
the mass splitting is large enough (5m ^> O(10) keV), the DM-nucleon inelastic scattering 
is suppressed kinematically. Thus, the EW-IMPs avoid the constraints from the direct 
detection experiments. On the other hand, large Sm may also induce sizable contribution 
of the Higgs-boson exchange to the DM-nucleon scattering at tree level. To avoid all these 
constraints, in this article, we simply assume O(10) keV <C 5m <C AM so that the elastic 
scattering of EW-IMP with nucleon is dominated by loop diagrams. 

In the following discussion, we consider the scenarios with Majorana fermion EW- 
IMPs, with either Y = or Y ^ 0. We simply express it as x° an d the charged Dirac 
components which couple with x° as i }± - 

It is straightforward to extend our calculation to scalar EW-IMPs. Contrary to the 
fermionic EW-IMPs, the scalar EW-IMPs have renormalizable self-interactions with Higgs 
boson, which contribute to the SI cross section of scalar EW-IMP with nucleon. When 
the interaction is suppressed enough to be ignored, elastic scattering of the EW-IMP with 
nucleon is induced by loop diagrams. In that case, the SI cross section of scalar EW-IMP 
with nucleon should agree with that of the fermionic EW-IMP when the EW-IMP mass is 
much larger than weak gauge boson masses. This is because the Lagrangians for both of 
fermionic and scalar EW-IMPs are the same in the non-relativistic limit or at the leading 
of the velocity expansion. Needless to say, the SD cross section of scalar EW-IMP with 
nucleon vanishes since the scalar EW-IMP has spin zero. 

3 Effective Interactions for DM-Nucleon Scattering 

In this section, we briefly review evaluation of the elastic scattering cross section of DM 
with nucleon, which is an important quantity for the direct detection experiments, as- 
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suming the DM is a Majorana fermion. The formulation given here is originally derived 
in Ref. [13]. See also Refs. [HI [T5] for further details. 

First, we write down the effective interactions of DM with light quarks (q = u, d, s) 
and gluorjf|: 

c cs = Yl c f + c f > ( 4 ) 

q=u,d,s 

Cf = d q X 7 M 75X° qivl5q + f q m q x°X° QQ 

+ + (5) 

Cf = fc¥x G%G a ^ . (6) 

Here, M and m q are the masses of the DM and quarks, respectively. The field strength 
tensor of the gluon field is denoted by G a . The last two terms in Eq. (JSJ) include the 
quark twist-2 operator, O q , which is defined as, 

^qi yDtfy + Dvjp - ^g^Pj q • (7) 

The coefficients of the operators in Eqs. (J5]) and (EJ) are to be determined in the following 
section. 

The cross section of DM with nucleon (N = p, n) is calculated from the effective 
Lagrangian in terms of SI and SD effective couplings, /jy and a N , 

a N = -m 2 R [\f N \ 2 + 3\a N \ 2 ] , (8) 

7T 

where m# = Mm^ j[M + m^) {m^ is nucleon mass). The SI and SD effective couplings 
are given as scattering matrix element of the effective Lagrangian between initial and final 
states. The results are 

f N /m N = £ f q f Tq + ^(g(2) + g(2))(^)+4 2 ))-^-/ TG / G , (9) 

q=u,d,s q=u,d,s,c,b 

£ d q Aq N , (10) 

q=u,d,s 



^ [IV 



a N 
where 



(N\m q qq\N)/m N = f Tq , f TG = 1 - £ f Tq , (11) 

u,d,s 

1 f 1 

m N 4 

(N\qwq\N) = 2s^q N . (13) 



(N(p)\Ol\N(p)) = —{p,Pu ~ -mWv) (q(2) + q(2)) , (12) 



3 We only keep the operators which give rise to the leading contributions for the DM-nucleon scattering 
with the non-relativistic velocity. Moreover, in order to remove the redundant terms, we use the integra- 
tion by parts and the classical equation of motion for the operators when we construct the low-energy 
effective Lagrangian [IB] . 
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Mass fraction 










(for proton) 




Second moment at \x = mz 


fau 


0.023 








(for 


3roton) 


fad 


0.032 




u(2) 


0.22 


u(2) 


0.034 


fas 


0.020 




d(2) 


0.11 


d(2) 


0.036 


(for 


neutron) 




s(2) 


0.026 


5(2) 


0.026 


fau 


0.017 




c(2) 


0.019 


c(2) 


0.019 


fad 


0.041 




6(2) 


0.012 


1(2) 


0.012 


fas 


0.020 











Spin fraction 


(for 


proton) 


Au p 


0.77 


Ad p 


-0.49 


As p 


-0.15 



Table 1: Parameters for quark and gluon matrix elements used in this paper. 

Here a s = gl/^n (g s is the coupling constant of SU(3)c), q(2) and q{2) are the second 
moments of the parton distribution functions (PDFs), and s M is the spin of the nucleon. 
Note the factor l/a s in front of fa in Eq. ([9]). It makes the gluon contribution comparable 
to the light-quark contribution, despite the interactions of DM with gluon are induced by 
higher- loop processes than those with light quarks j8]. In the present CclSC, clS discussed 
in the previous section, the DM-quark tree-level scattering is highly suppressed, and thus 
the one-loop processes become dominant for the DM-light quark effective interactions. 
Therefore, in order to give accurate prediction of the scattering cross section, we should 
evaluate not only the one-loop diagrams with light quarks but also the two-loop diagrams 
with gluon. 

In Table 1 we list the numerical values for the parameters that we used in this article. 
The mass fractions of light quarks, fa q (q = u,d,s), are calculated using the results in 
Refs. [I7J dSl [T9] . The procedure for evaluating them is described in these references and 
Ref. [14J. The second moments of PDFs of quarks and anti-quarks are calculated at the 
scale /i = mz (mz is Z-boson mass) using the CTEQ parton distribution [2oj^|. Finally, 
the spin fractions, Aq^, in Eq. (fT3"l) are obtained from Ref. [22]. The second moments and 
the spin fractions for neutron are to be obtained by exchanging the values of up quark 
for those of down quark. 



4 Results 

4.1 Coefficients of the Effective Lagrangian 

Now we evaluate the coefficients of the effective interactions displayed in Eqs. ([5]) and (J6j). 
From Eqs. ([2]) and (J3j) , it is found that the EW-IMP, x°, interacts with the weak gauge 

4 As will be described later, the terms with quark twist-2 operators in Eq. ([5]) are induced by the 
one-loop diagrams in which the loop momentum around the weak-boson mass scale yields dominant 
contribution. This fact verifies the use of the second moments of PDFs at the mz scale. See also 
discussion relevant to it in Ref. [21] . 
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Figure 1: One- loop diagrams which induce effective interactions of EW-IMP DM with 
light quarks. There are also W-(Z-) boson crossing diagrams, which are not shown here. 




Figure 2: One- loop diagrams which correspond to the one-loop quantum correction to the 
EW-IMP-Z (7) interaction vertex. These contributions turn out to vanish. 



bosons (W*, Z®) as 

192 ry> , nv , 1X , nr+ , 92 



AC 



hit . 



jV 7 ™ 2 - ( 2Y + !) 2 X°7^" + jv/n 2 - (2Y- l) 2 xVV> + W" + h.c. 



COStV 



(14) 



where is the weak mixing angle. The Majorana field 77° is introduced for the cases of 
7^0. (See Eq. ©.) In either case (Y = or Y ^ 0), the EW-IMP does not have any 
interaction by itself. Thus, it is loop diagrams that yield the leading contribution to the 
EW-IMP-nucleon elastic scattering cross section. 

First, we consider the one- loop processes. The relevant diagrams are shown in Figs. [1] 
and 121 The diagrams in Fig. [T] give rise to the coefficients in Eq. (J5J) as 
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4m 2 



n 2 - (AY 2 + 1) 
8m w 



gn(,w) + 



Y 2 



Amzcos^Ow 



9h{z) 



+ 



(«) 2 - «) 2 ) ^ 2 a£ 



cos 4 # 



w 



rrrr 



■gs{z), 



d„ 



4 2) 



n 2 - (AY 2 + 1) a 2 , . , 2 («) 2 + «) 2 ) E 2 



m 



— 0av(w) + 



(4E 2 + i) « 2 



m 



■gm(w) + 



n 



w 

2 



(AY 2 + 1) a 2 

-— f-^H + 



m 



COS 4 ^^/ 


2 i 

774 


2 («) 2 + « 


) 2 ) E 2 


a 2 


COS 4 #^ 






2 («) 2 + « 


) 2 ) Y 2 


a 2 


COS 4 ^ 




■? 



9T2{Z). 



T 3q - Q q sin 2 6> 



w ■ 



-ir. 



3q • 



(15) 
(16) 
(17) 
(18) 



Here, and are the masses of Higgs and W bosons, respectively, and a 2 = g^/Air. 
We also define the vector and axial-vector couplings of quarks with Z boson as 



(19) 



where T 3q and Q q denote the weak isospin and the charge of quark q, respectively. Fur- 
thermore, we parametrize w = m^/M 2 and z = m|/M 2 in the above expressions. The 



first term in Eq. (1151) is induced by the Higgs-boson exchange process, shown in diagram 
(a) of Fig. [TJ The other terms in Eqs. ( fl5lfT8l) are all obtained from diagram (b). The 
mass functions, Qh(?z), 9s{ x )i 9av(x), 9ti(x), and gr2(x), in Eqs. f fT5lfT8|) are given in Ap- 
pendix |A] We ignored the mass differences between x° an d ip and also between x° and 
f}° here. The loop integrals for the diagrams are finite, and they are dominated by the 
loop momentum around the weak-boson mass scale. We note here that some of these mass 
functions do not vanish in the limit of w, z — >■ [5]. This implies that the SI interactions 
of EW-IMP with light quarks are not suppressed even if the EW-IMP mass is much larger 
than the gauge boson masses as we described in the Introduction. On the other hand, we 
found that both diagrams in Fig. [2] vanish separately by explicit calculation. 

Next, let us move to the two- loop diagrams which yield the effective scalar coupling of 
EW-IMP with gluon, fa- These diagrams are presented in Fig. [3J For each diagram, the 
gluon contributions are classified into two types in terms of the momentum scale which 
dominates in the loop integration. So-called "short-distance" contribution means that the 
momentum is typically masses of heavy particles, such as DM particle or the weak gauge 
bosons in the diagrams, and "long-distance" one means that the momentum is typically 
mass of quark in the loop diagrams. Among the latter one, the contributions in which the 
light quarks run are already incorporated in the mass fractions defined in Eq. ffTTj) . 
Therefore, we do not need to add them in the calculation of the gluon contribution; 
otherwise we would doubly count them [TJ]. Consequently, the gluon contribution from 
each diagram is given by 



f (0 

JG 



E 



f(*)|SD i 
JG \q 



=u,d,s,c,b,t 



E 

Q=c,b,t 



LD 

Q ■ 



(20) 
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Figure 3: Relevant two-loop diagrams which contribute to effective scalar coupling of 
EW-IMP DM with gluon. There are also W-(Z-) boson crossing diagrams, which are not 
shown here. 



Here /^l, and /^|g D denote the short-distance and long-distance contributions of 
quark q in the loop in diagram (i) (i = a, b, c) of Fig. [3j respectively. We also take 
into account large QCD corrections in the long-distance contributions [23] by using cq = 
\ + \\a s {mQ) / An (Q = c, b, t). We take c c = 1.32, q, = 1.19, and c t = 1 for a s {mz) = 0.118 
in our calculation. Note that the long-distance contribution is gauge invariant. This is 
because its contribution to the operator faG a G a ^ u is evaluated from scalar-type effec- 
tive operator f q m q qq jT3j [H]. (See also later discussion where the explicit calculations 
are given.) Consequently, the gauge invariance of the short-distance contribution is guar- 
anteed since summation of the both contribution is obviously gauge invariant. Then, the 
effective coupling of EW-IMP with gluon is obtained as 

fa = f ( G ] + f { G ] + f ( G ] ■ (21) 

Let us see each diagram closely. It is obvious that diagram (a) gives the long-distance 
contribution. Thus, we sum up for heavy quarks in the loop, and get 



,(o) _ a s «2 \^ 

/g " 12tt X Am\ ^ ° Q 

h Q=c,b,t 



2 (A\r2 

(22) 



n„ n- \ ^ ~n 2 - {AY 2 + 1) Y 



9h(w) + — — g H (z) 

mw AmzcosWw 



Here the first and second terms in the bracket come from W- and Z-boson exchanges, 
respectively. As we described above, this long-distance contribution is given by effective 
scalar-type coupling (of the Higgs contribution) as —j^fo.- 

For the calculation of diagrams (b) and (c), on the other hand, we follow the steps 
supplied in Ref. [13] . In the work, the systematic calculation for the VT-boson exchange 
diagrams at two-loop level in the Wino DM scenario is given. The procedure is applicable 
to compute the two-loop diagrams in our case. For the H^-boson exchange diagrams, the 
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analytic result is simply given by multiplying a factor (n 2 — (AY 2 + l))/8 to their result 
of the Wino DM case, (as the first term in Eq. (122j) ). Even for the Z-boson exchange 
contribution, the calculation is straightforward extension. In the following we outline the 
calculation and give the result of the Z-boson exchange contribution. 

First, we calculate the vacuum polarization tensor of Z boson in the gluon background 
field, n^tjo). It is given as 



KM 



= - E / (2 t)w% TrL + c W< + <^)s q {i)iM v q + <75 W - p)h (23) 

where TrL+c denotes the trace over the Lorentz and color indices, and [q] means the 
sum is taken over all quarks for the short-distance contribution and heavy quarks for the 
long-distance contribution. The quark propagators S q (p) and S q (p) are under the gluon 
background field with the Fock-Schwinger gauge. (See Appendix A in Ref. [14] .) We 
decompose the polarization function as 

KM = (-^ + Pj f) n|(P 2 ) + v -f- nf(p 2 ). (24) 

By the explicit calculation, we found that the longitudinal part IIlQo 2 ) does not contribute 
to /(Jj. Thus, we only calculate transverse component, IT^(p 2 ). The contributions from 
diagrams (b) and (c) to Ilf (p 2 ) are given as 



-2(a*) 2 (p 2 B^ + 4Bir%(25) 



where [TlfQo 2 )]^ denotes the contribution from diagram (i) [i = b,c) with quark q 
running in the loop. Here, u \gg" represents that we keep only the terms which are 
proportional to G^ u G afJ-l/ in the polarization function, since the terms proportional to 
G° G"'" contribute to fa- The loop functions in Eq. 025]) are defined as follows: 

H B (n,m) = f d_k k» 

1 1 J m 2 [k 2 - m 2 q ] n [{k + p) 2 - m 2 q ] m ' 1 J 

ppB 21 +g JJ 22 - J m2[k2 _ m2]n[{k + p)2 _ m2]m ■ V<) 



5 This is consistent with the fact that EW-IMPs have no tree-level coupling with Higgs boson and 
Goldstone bosons of the electroweak symmetry breaking, which turn into the longitudinal modes of Z 
and W bosons. 
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As is the same in the calculation of two-loop W-boson exchange diagrams discussed in 
Ref. [H], the diagrams (b) and (c) yield the short-distance contribution and long-distance 
contribution, respectively. We have checked this identification by explicit calculation. 
This is also confirmed from the fact that the long distance-contribution corresponds to the 
one which is calculated from quark triangle diagram using effective scalar-type coupling 
f q in the M and mw/z oo limit [14], as we mentioned before. Therefore, ILp (p 2 )\gg is 
obtained as 



nf(p 2 )| GG = £ [nf( P 2 )]; ) | GG+ £ CQ [n|( P 2 )]«| GG 

q=u,d,s,c,b,t Q=c,b,t 



(28) 



The results of [flf (p 2 )]^ \ GG and [U§(p 2 )]" c) | 

gg are given in Appendix 



With these Z-boson polarization functions, we get the contribution from Z-boson 
exchange diagrams at two-loop level. Combining it with the contribution from VF-boson 
exchange diagrams, we eventually obtain the effective coupling f G from diagrams (b) and 
(c) as 



(<0 



a s at 

47T 



n 



{ay 2 + r 



9w{w,y) 



Y' 2 



w 



4m|cos 4 6 l vK 



9z(z,y) 



(29) 



where y = m 2 /M 2 {m t is the top quark mass), and the mass functions, gw(z,y) and 
gz(z, y), are given in Appendix |B] We ignore the mass of quarks except that of top quark 
in this computation. 



4.2 Spin-independent and Spin-dependent Cross Sections 

With the effective couplings derived above, we evaluate the cross section of EW-IMP 
with nucleon. First, we discuss the SI cross section. In order to look over the behavior 
of the cross section, we plot the SI cross section as a function of the EW-IMP mass in 
Fig. HI In this figure, we set n = 5 and Y = 2 as an example of the case where both 
W- and Z-boson exchange diagrams contribute, and take Higgs-boson mass as =130, 
115, 300 and 500 GeV from bottom to top. It is found that the SI cross section has little 
dependence on the EW-IMP mass when the mass is larger than 0(1) TeV. We have 
checked such a behavior for other n and Y cases. The cross section ranges from 10~ 47 cm 2 
to 10~ 46 cm 2 in this figure when Higgs-boson mass increased from = 115 GeV to 
500 GeV. (Although it is not monotonic.) 

In Fig. we plot the EW-IMP-proton SI cross section as a function of Higgs-boson 
mass for n = 5 (left panel) and both n = 2 and n = 3 (right panel). The EW-IMP mass 
is taken to be equal to 1, 2.7, and 10 TeV for n — 2, 3, and 5, respectively^ [HI HP] - It is 
seen that the cross section is enhanced as n is larger and Y is smaller. This behavior will 

6 Recall that the EW-IMP mass M — 1, 2.7, and 10 TeV is suggested in order to explain the observed 
DM abundance when n — 2, 3, and 5, respectively, and Y = 0. The thermal relic abundance of the 
EW-IMP has not been evaluated when Y ^ 0. We use those suggested values in Y ^ case since the SI 
cross section is insensitive to the EW-IMP mass as mentioned in the text. 
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J I I I I I I I I I I I 1—1— 



1 10 
DM mass [TeV] 

Figure 4: DM-proton SI cross section for n = 5 and Y — 2. We take Higgs-boson mass 
as rrih =130, 115, 300 and 500 GeV from bottom to top in this figure. 



be explained later. We also found that the SI cross section of EW-IMP with nucleon is 
far below the current experimental bound. This is the consequence of the calculation in 
which all the relevant terms at leading order are taken into account. The suppression of 
the cross section originates in an accidental cancellation within the SI effective coupling, 
/at. Such an cancellation was already pointed out for the Wino dark matter, i.e. n = 3 
with Y = [8j. In our calculation, we found that the similar cancellation also occurs in 
general n-tuplet cases with both W^-boson and Z-boson contributions. 

Let us examine what happens in the SI effective coupling. In Fig. [6j we show the 
contribution to f p /m p of each effective coupling in Eq. The left (right) panel in 

the figure shows each contribution due to the W-boson (Z-boson) loops as a function of 
Higgs-boson mass. The plot is given in units of (n 2 - (AY 2 + l))/8 and Y 2 in left and 
right panels, respectively. In the figure, we show the contribution of the Higgs-boson 
exchange including both one- and two- loop contributions (solid), the gluon contribution 
except for the Higgs-boson contribution (dashed), the quark twist-2 operator contribution 
(dash-dotted), and the contribution to the quark scalar-type operators with the coefficient 
f q from the Z-boson box diagrams (double-dashed). The VT-boson box diagrams do not 
generate the scalar-type operators [8]. It is found that the contributions from the twist-2 
operators generated by both VT-boson and Z-boson loops are positive, while all the other 
contributions are negative. In addition, they have roughly the same order in absolute 
value. This is why the cancellation happens in the SI effective coupling induced by 
each iy-boson and Z-boson contribution. We also find that the Z-boson contribution in 
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Figure 5: SI cross sections of DM-proton elastic scattering as a function of Higgs-boson 
mass for n = 5 (left panel) and n = 2, 3 (right panel). We take M = 1, 2.7, 10 TeV for 
n =2, 3, and 5, respectively. In left panel, solid, dashed, and dash-dotted lines represent 
n = 5 with Y — 0, 1, and 2 cases, respectively. In right panel (n, Y) = (3,0), (3, 1), and 
(2, 1/2) correspond to solid, dashed, and dash-dotted lines, respectively. 
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Figure 6: The contributions of each effective coupling in Eq. (Q to f p /m p . Left panel 
shows each contribution from P^-boson loops as a function of Higgs-boson mass, while 
right one illustrates contribution by Z-boson loops. These plots are shown in units of 
(n 2 - (AY 2 + l))/8 and Y 2 , respectively. We take M = 10 TeV. In the both panels, 
solid line represents the contribution of the Higgs-boson exchange (including both one- 
and two- loop contribution), dashed line indicates the rest of gluon contribution in the 
two-loop contribution, and dash-dotted line denotes that from quark twist-2 operators. 
Double-dashed line in right graph shows contribution of scalar-type operator coming from 
box diagrams. 
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Figure 7: SD cross section of DM-proton elastic scattering for n = 5 (left) and n = 2, 3 
(right). In left graph, solid line, dashed line, and dash-dotted line represent Y = 0, 1, 2, 
while in right graph they represents (n, Y) = (3,0), (3,1), (2, 1/2), respectively. 

units of Y 2 is generally rather small compared with the W^-boson contribution in units 
of (n 2 - (AY 2 + l))/8. Consequently, the cross section is suppressed when Y gets larger 
(because the coupling between W boson and EW-IMP is suppressed). When Higgs-boson 
mass is relatively small as rrih ~ 115 — 200 GeV, the contribution of Higgs-boson exchange 
becomes large. Then, there is much more significant cancellation within each W- and Z- 
boson contribution. As a result, these cancellations give rise to a very suppressed SI cross 
section. 

Next, we show the SD cross section in Fig. [7] for completeness. Since the SD cross 
section is independent of Higgs-boson mass, we plot it as a function of the EW-IMP mass. 
Contrary to the SI cross section, the SD cross section is suppressed when the EW-IMP 
mass is larger [5]. In the left panel, the SD cross sections of n = 5 cases with Y — 0, 1, 
and 2 are shown in the solid line, the dashed line, and the dash-dotted line, respectively. 
In the right panel, we give the cross sections by taking (n, Y) = (3, 0), (3,F) and (2, 1/2), 
in the solid line, the dashed line, and the dash-dotted line, respectively. \j As expected, 
the cross section is enhanced by larger n in a similar way to the SI cross section. The 
figure also shows that the cross section decreases when one sets Y larger, because the 
contribution from W^-boson exchange reduces. We also find that the contribution from 
Z-boson exchange has the opposite sign and rather small significance compared to that 
from H^-boson loop diagrams. 

Throughout the calculation, we have ignored the mass splitting 8m between the neutral 
components for case. The magnitude of 8m is model-dependent. If 8m is relatively 

large against our assumption, the Higgs-boson exchange contribution at tree level becomes 
significant, and one might expect a larger SI cross section. Even in such a case, the 

7 There was an numerical error in the result of SD cross section in Ref. [8] . We have corrected it here. 
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extension is straightforward. 

Lastly, we refer to the difference between our results and previous works. It is found 
that the mass functions in Eq. (|30|) in the Appendix [A] which gives scalar- and twist- 
2-type (and also axial-type) contributions agree with the corresponding loop functions 
given in Eq. (42) of Ref. [5], except for F^^x) and F^(x) in Ref. [5]. (This was pre- 
viously described in Ref. [8j. ) On the other hand, gluon contribution was taken into 
account without explicit calculation; however, the sign of the contribution is the same 
with our result. As a consequence, similarly to our result, the cancellations between each 
contribution happen and the cross sections computed in their works are comparable to 
our results. 

On the other hand, the SI cross sections in Refs. [61 [7J are larger than our results by 
more than an order of magnitude. In Ref. [5] only the scalar-type operator is taken into 
account in the limit M 3> my/ 3> m q , which does not agree with our result in the same 
limit. Although they mention the twist-2 type contribution, they took into account this 
contribution in the same sign as the scalar-type one. Moreover, the gluon contribution is 
omitted. Thus, there is no cancellation which we found in the effective coupling, to give 
rise to the large cross section. 

In Ref. [7j, the explicit results of the loop calculation are given. However, all of the 
loop functions except ff(x) in Eq. (26) in Ref. [7J are different from ours. To check their 
result further, we compare them with ours in the limit mw/z/M — > 0. In Eq. ( TT5|) . which 
is scalar-type contribution, the first term (of the first parenthesis) agree with their result 
except for the sign, and the other terms are —1/2 of corresponding term in their work. 
For twist-2 type contributions, the one (i.e., Eq. (1T7|) ) is the same as their result in the 
limit, while another contribution (i.e., Eq. (ITS]) ) is neglected in their paper. The gluon 
contribution is considered, but not evaluated properly. As a result, the cancellation we 
found in our calculation does not occur, which makes the resultant cross section large. 

5 Conclusion 

We have studied the direct detection of EW-IMP DM, which is the lightest neutral com- 
ponent of an SU(2)l multiplet interacting with quarks and leptons only through the 
SU(2) L x U(1)y gauge interactions. Although EW-IMP DM does not scatter off nucleon 
at tree level, it does at loop level, and the SI scattering cross section is not suppressed 
even if the mass of EW-IMP is much larger than gauge boson mass. We evaluate the 
two-loop processes for the EW-IMP-gluon interaction in addition to the one-loop pro- 
cesses of EW-IMP scattering with light quarks, since both of them yield considerable 
contribution to the SI scattering cross section. As a result, the SI cross section is found 
to be 0(1O~^ 46_48 )) cm 2 , depending on Higgs-boson mass, the number of components n in 
the multiplet, and the hypercharge Y of the EW-IMP DM. Such small value of the cross 
section is due to an accidental cancellation in the SI effective coupling. This cancellation 
is a general feature for EW-IMP DM, and thus makes it difficult to catch EW-IMP DM 
in the direct detection experiments in the near future. 
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Appendix 

In this Appendix we provide the mass functions and the Z-boson vacuum polarization 
functions, both presented in Section HI 

A Mass Functions from One-loop Diagrams 

The mass functions used in Eqs. ( jl~5]fT8]) in the calculation of one-loop diagrams are 
gn(x) = -^(2 + 2x-x 2 )tan' 1 ] + 2^(2 - x log(x)) , 

Ox \V X J 

gs(x) = ^(4-2x + a; 2 )tan- 1 f^j + ^(2 - x log(x)) , 
9av(x) = ^-^(8 - x - x 2 ) tan -1 f^=J ~ ^ x ( 2 ~ (3 + a;) log(x)) , 
9ti(x) = ^6,(2 + x 2 ) tan^ 1 ^^^j + ^^(1 - 2x - x(2 - x) log(x)) , 
g T2 (x) = -^j-x(2 - 4x + x 2 ) tan^ 1 (rj=J - ^Vx(l - 2x - x(2 - x)\og(x)) , (30) 
with b x = a/1 — x/4. 

B Mass functions from Two-loop Diagrams 

Here, we present the polarization functions and the mass functions for Z boson. Perform- 
ing the integration in Eq. (l25]) . we derive the polarization functions as 

2(p 2 - 4mJ)(p2 - 3ml) + im W " em*)^ - lcotr 1 Lffl- - l) 

p 2 (p 2 — Am 2 ) 2 
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«) 2 (2p 4 - limy + 24m*) - (af) 2 {p 2 - Am 2 q ){2p 2 ~ 3m 2 q ) 



p 2 (p 2 — Am 2 ) 2 



p 2 {p 2 — Am 2 ) 3 



and in the limit of zero quark mass, these functions lead to 



[n?b 2 )f 



(6) 



|GG 



3 cos 2 6 1 



:(G^) 2 [( 



V\2 /„A\2 



(31) 



(32) 



Using the polarization functions, we compute gz(z,y) in Eq. fT2Ul) as 



9z{z,y) 



£ {«) 2 + «) 2 }- 2 £ c q{( 

,q=u,d,s,c,b Q=c,b 



a V ) 2_ (a A ) 2* 



x 4£bi(z) +g t (z,y) . 

(33) 

Here, the first term comes from the contribution of all quarks except top quark, and cq 
in it represents the QCD corrections in the long-distance contributions [23]. The second 
term g t (z,y) is the contribution of top quark. The function <?bi(;e) is given as 



1 (x 2 -2x + 4)tan- 1 ( 2 ^) 
9bi(x) = -— y/x(x\og(x) - 2) + — ^- 



(34) 



which is equal to the one in Ref. [8]. In the calculation of g t (z,y), we decompose it into 
two parts: 



9t(z, y) = 9t°' log (z, y) + g? s (z, y). 



gt°~ og (z,y) is analytically obtained as 



where 
Gn(z,y) 



9 r l ° s (z,y) = (aTf G n (z,y) + (a?) 2 G t2 (z,y), 
sfz(Yly 2 - zy + z 2 ) 



(35) 



(36) 



+ 



3(4y - zf 
z 3/2 (A8y 3 - 20zy 2 + Ylz 2 y - z 3 ) 

6(4y-zf ' 3(4y-zf 

z 3 / 2 ^/y(lQy 3 - 4(2 + 7z)y 2 + 14(2 + z)y + 5z) _ x ( ' VT=y 



2z 3 l 2 y 2 Uy-lz) , lA s 
log ^ + J. K y NO ; log(4y) 



3(4y - zf^T^y- 



tan 



tan 



48(z 2 -2z + A)y 3 - 4z(5z 2 - 10 z + U)y 2 + \2z 3 {z - 2)y - z 3 (z 2 - 2z + 4) 

3(4?/ - z) 3 y/£=z : 
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y/i(2y - z) z^jSy'-Szy + z 2 ) Az 3 / 2 y 2 , (A , 

Gt2{z ' y) = ~W^) — logz ~ W^f log(%) 

4z 3 / 2 yy(2y 2 - y - 1) _ t /VT=y\ 

+ (i^vr^ tan l^rj ^ 

8z(z 2 - 2z + l)y - (z 2 - 22 + 4)(8y 2 + z 2 ) _ x ( 



(% - z)V4 - i V v 7 ^ 

On the other hand, we calculate g\° s (z,y) numerically. For convenience, we rewrite it as 

g^(z, y) = Az^ 2 y 2 (A, y [h + I 2 ] + A 2 [J a + h]) , (38) 

with 



A 1 = -2(ar) 2 + 4(af) 2 , 
A 2 = -(af) 2 + (af ) 2 , 

and then carry out the following integrals numerically, 

J = r° dt + y^) ( log [y/t + ^ + y/t] - log [y/t + %- v^]) 



(39) 







[t + z] 2 [t + Ayf/ 2 t 

h= i°° dt l x (t + 2-VtVt + 4) H [Vt + 4y-+Vt\ -log[vT+4^- V*]) 
/ 3 = / dt 







2 [t + z} 2 [t + 4y] 5 / 2 *V2 

(y/t + 4 - Vf) (log [yT+% + Vt\ - log [y^T% - Vj]) 
[t + z] 2 [t + 4yf/ 2 



J4 — / dT x 







2 [t + z] 2 [t + %] 5 /2 



Lastly, we present the mass function g^(w,y) in Eq. ( 1291) . which is readily obtained 
by following the similar procedure described in Ref. [8]: 

9w(w, y) = 2g B1 (w) + g B 3(w, y) (40) 

with 

W x > v) = 9b3( x > y) + Cb9ial( x > y)- ( 41 ) 

Here the first term of g^(w,y) is coming from the first- and second-generation quark 
loop diagrams, while the second term is from the third-generation quark loop diagrams. 
The mass function gBi(x) is displayed in Eq. ( 1341) . Although we use the same symbol for 
gBs(x,y) as in Ref. [8], it is reevaluated with the QCD correction for the long-distance 
contributions, which we illustrate with the factor q, in Eq. ( |4~TT) explicitly. We have checked 
numerically that including the QCD correction changes the SI cross section by up to a 
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few %. The functions g^l(x } y) and g^(x,y) are analytically given as 

(i) -x 3 / 2 -x 3/2 y 2 x 5/2 (x-2y) 

9 ™ (X ' V) = Y2^x) + 2^y—^ hgV ~ 2A(y-x) 2 hg 
x 3 / 2 ^y(y + 2)^A^ _ 1 f ^=j[ 

: r~r tan I 

io/», „\2 \ nr. 



12(2/ -xf V Vv J 

x(x 3 -2(y+ l)x 2 + A(y+ l)x+4y) _ x 

~\~ — tclll 

12(y - x) 2 V4 - x 

(2)/ v _ -x 3/2 y -x b/2 y 2 x 5/2 y 2 . 

9B3 [x, V) ~ 12{y _ x? + 24{y _ x)3 1o § y + 24 ( y - x) 3 



/x 



(42) 



respectively. 



12(7/ - + 24(y - x) 3 l0g2/ + 2A(y - a;) 3 ^ 
x 3/2 y/y(-6y + ^ 2 - 2xy - 2x) _ 1 A/3 - y > 
12(y-a;)3V^ ^ V > 

—xy{x 2 y — 2xy — 6x — 2j/) _ x / \/4 — s\ 
12(j/-x) 3 v 7 ^"? an V v 7 ^ / ' 
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